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We study the energetics, the electronic structure, and the ballistic transport of an infinite Au 
monatomic chain with an adsorbed CO molecule. We find that the bridge adsorption site is ener- 
getically favoured with respect to the atop site, both at the equilibrium Au-Au spacing of the chain 
and at larger spacings. Instead, a substitutional configuration requires a very elongated Au-Au 
bond, well above the rupture distance of the pristine Au chain. The electronic structure properties 
can be described by the Blyholder model, which involves the formation of bonding/antibonding 
pairs of 5cr and 27r* states through the hybridization between molecular levels of CO and metallic 
states of the chain. In the atop geometry, we find an almost vanishing conductance due to the 5a 
antibonding states giving rise to a Fano-like destructive interference close to the Fermi energy. In 
the bridge geometry, instead, the same states are shifted to higher energies and the conductance 
reduction with respect to pristine Au chain is much smaller. We also examine the effects of strain on 
the ballistic transport, finding opposite behaviors for the atop and bridge conductances. Only the 
bridge geometry shows a strain dependence compatible with the experimental conductance traces. 

PACS numbers: 73.63.Rt, 73. 23. Ad, 73.20.Hb 



C^ 



o 

> 

rn 
O 
(N 



X 



I. INTRODUCTION 

Au monatomic chains are produced in nanocontact 
experiments through different techniques ji such as the 
Scanning Tunneling Microscope,—!^ the High-Resolution 
Transmission Electron Microscope)^ and the Mechani- 
cally Controllable Break- Junction (MCBJ).^ It has been 
demonstrated that a single row of suspended Au atoms 
between two electrical contacts typically displays a con- 
ductance of one quantum unit Gq = 2 e^ jh (e being 
the electron charge and h Planck's constant),?^ which can 
be satisfactorily explained by a ballistic transport model 
based on the Landauer-Biittiker theory and a single, al- 
most perfectly transmitted conductance channel (as also 
confirmed by shot noise experiments, see Ref. |l|)- This 
conductance channel arises from the s-type valence elec- 
tronic states of the Au chain, which form a wide, spin- 
degenerate electronic band crossing the Fermi level. 

The reduced dimensionality of monatomic chains, and 
of atomic-sized contacts in general, favors the adsorp- 
tion of impurities on the atoms forming the thinnest 
part of the contact because of their reduced coordination 
with respect to bulk or surface atoms i^"— It is nowadays 
well understood that accidental or controlled contami- 
nations with small impurities may change substantially 
the nanocontact conductance, because an adsorbed im- 
purity can act as a source of electron scattering in the 
contact region. On the experimental side, a large num- 
ber of cases have been already examined so far, such 
as H2 on Ptji£"— Auji^iii and Pd nanocontactsr^ but 
also on Fe, Co, and Ni nanocontacts;^^ CO on Pt,— i^i 
Au, Cu, and Ni nanocontacts^^ O2 on Au and Ag 
nanocontactS)^ and many others. The effect of contam- 



ination is usually visible in the conductance histograms: 
the distribution of peaks often looks significantly altered 
after the admission of impurities in the proximity of the 
nanocontactjiSiHii^iiSrJ^ In some cases, the presence of 
adsorbed impurities changes the formation process and 
the structural properties of the chain: a remarkable ex- 
ample is the enhancement of monatomic chains formation 
for Au nanocontacts in 02-enriched atmosphere, which 
has been attributed to atomic oxygen embedded in the 
chains J^iiS, 

Several theoretical studies have been carried out in this 
field,— '^"^1^^" — but the understanding of the physical prop- 
erties of monatomic chains in presence of impurities is 
far from complete. In this respect, an interesting case 
study is the adsorption of CO on Au monatomic chains. 
Recently, it has been demonstrated experimentally that 
CO can interact with monatomic chains produced in Au 
nanocontacts, modifying their conductance.^^ However, 
compared to transition metal nanocontacts (Pt and Ni, 
for instance), CO only marginally alters the conductance 
histograms of Au nanocontacts. CO adsorption on Au 
chains has also been demonstrated for short monatomic 
chains grown on a NiAl substrate;^ but with a much 
stronger effect on the electronic states of the Au chain. 
Indeed, both the experimental results^ and the elec- 
tronic structure calculations^! indicate that the delocal- 
ization of some electronic states of the chain can be de- 
stroyed by a CO molecule adsorbed in the atop posi- 
tion. According to ballistic conductance calculations for 
an unsupported chain ;^^'^^ the impurity in the atop posi- 
tion acts as a "chemical scissor" for the coherent electron 
transport along the chain, resulting in a strong depression 
of the electron transmission near the Fermi level. These 



theoretical results are clearly not adequate to satisfacto- 
rily explain the interaction of CO with Au monatomic 
wires and Au nanocontacts, since the conductance his- 
togram of Au in presence of CO does not display any low 
conductance tail.^^ The large strain imposed by the NiAl 
substrate on the Au chain, corresponding to a Au-Au 
spacing close to the rupture point of the isolated chain, 
and the atop position of CO might not be representative 
of the real adsorption geometry on suspended chains in 
Au nanocontacts. Indeed, another theoretical study sug- 
gests that an adsorption geometry with a CO molecule at 
the bridge site of a short Au chain between surfaces gives 
only a very small reduction of the ballistic conductance 
with respect to the clean Au nanocontact.— 

In this work, we carry out a density functional study of 
the interaction between CO and a monatomic Au chain, 
focusing on the bonding mechanism and its effects on 
the electron transport properties of the chain. The ad- 
sorption energetics of a single CO molecule at the bridge 
site or at the atop site are studied as a function of the 
Au strain. We also consider a linear geometry with a 
substitutional CO, because it has been suggested as a 
possible intermediate meta-stable configuration respon- 
sible for the conductance increase observed while pulling 
a Au nanocontact contaminated with CO."'^^ A simplified 
model is adopted here, consisting of an infinite, straight 
monatomic chain with a homogeneous spacing between 
the Au atoms (see Fig.[T]). We find that the bridge geom- 
etry is energetically favored with respect to the atop ge- 
ometry for all Au strains considered here, while the linear 
substitutional configuration could occur only in presence 
of a very elongated Au-Au bond. 

We show that the electronic structure of both the 
bridge and atop adsorption geometries can be inter- 
preted through the Blyholder modeP*^ (and its successive 
refinements'^^'^^), which is commonly used to explain the 
interaction of CO with transition metal surfaces and has 
already been proven valid in the case of CO adsorbed on a 
Pt monatomic chaiui^i^i We also find a significant corre- 
lation between the electronic transmission across the Au 
chain and the electronic structure, similarly to what seen 
for Pt chains in presence of COi^ For instance, a large 
reduction of the s-electron transmission is observed in 
correspondence of the density-of-states maximum of the 
5(Ja states. The energy and broadening of these states 
strongly depend on the adsorption site and on the Au 
strain, making the s transmission rather sensible to the 
atomic geometry. 

We also investigate how the ballistic conductance de- 
pends on the Au strain, finding two very different be- 
haviours in the bridge and atop geometries: in the for- 
mer, the conductance slightly increases with strain, while 
in the latter it progressively vanishes. This is explained 
by the different energy of the 5aa states, much closer to 
the Fermi level (Ep) in the atop geometry than in the 
bridge one, and by their energy shift with strain, toward 
E-p in the atop geometry but away from Ep in the bridge 
one. The present study is mainly aimed at a qualitative 
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FIG. 1. (Color online) Sideview of the atomic geometries used 
to model CO adsorption on the infinite Au chain: at the (a) 
bridge and (b) atop sites, and in the (c) linear substitutional 
and (d) tilted bridge configurations. The tetragonal supercell 
is indicated for the bridge geometry with dashed lines. The 
unit ceU for the metaUic leads (one Au atom) is shown on the 
left, enclosed by dashed lines. 



analysis of the ballistic conductance and will be succes- 
sively complemented by other work to address two points 
which are essential to give a more quantitative picture. 
First, the theoretical conductance at low Au strains suf- 
fers from a spurious contribution from d electrons due 
to their self-interaction error, which should be tackled 
within a beyond-LDA approach. ^^ Second, the calculated 
values of the tipless conductances of the bridge geometry 
are somewhat too large compared to the experimental 
values and thus call for a more realistic modeling of the 
nanocontact.— 

The paper is organized as follows: after describing 
in Sec. |TT] the numerical methods and approximations 
adopted in our calculations, in Sec. IIIII we shall present 
the energetics of CO adsorption on the infinite Au chain 
at different strains. In Sec. lIVI and Sec.|V]we will present 
the electronic structure and the strain-dependent elec- 
tronic transport properties, respectively, for the tipless 
Au chain in presence of CO. Finally, our conclusions will 
follow in Sec. IVTl 



II. 



METHODS AND COMPUTATIONAL 
DETAILS 



We carried out first principles calculations in the 
framework of density functional theoryS^ (DFT) solv- 
ing the single-particle Kohn-Sham (KS) equations^^ for 
a system of atoms in periodic boundary conditions as 
implemented in the electronic structure code PWscf , in- 
cluded in the Quantum ESPRESSO packagei^ We 
used the Perdew and Zunger*'^ parametrization for the 



exchange and correlation energy within the local den- 
sity approximation (LDA) and, additionally, we com- 
puted optimal distances and chemisorption energies in 
the generalized gradient approximation (GGA) using the 
Perdew-Burke-Ernzerhof (PBE) functional."*^ Ultrasoft 
pseudopotentials^ (USPPs) have been used to describe 
nuclei and core electrons. For Au, within the LDA we use 
the USPP described in Ref. 43, while within the GGA we 
generated a new USPP using the reference configuration 
and cutoff radii of the LDA USPP. For C and O, we use 
the same USPPs as in Ref.[3j. The KS wavefunctions and 
the charge density are described using a plane wave basis 
set with kinetic energy cut-offs of 32 Ry and 320 Ry, re- 
spectively. Electronic occupations are broadened using a 
standard smearing technique,-^ with a smearing param- 
eter a = 0.01 Ry. We partially relaxed our structures 
(see later in the text) with the BroydenFletcherGoldfarb- 
Shanno algorithm'^^ optimizing selected atomic positions 
until atomic forces get lower than 0.026 eV/A. 

To simulate the infinite chain model with one CO 
molecule we use the supercell method, as exemplified in 
Figs. [TJa) and [TJd) for the bridge, atop, substitutional 
and tilted-bridge geometries, respectively. Tetragonal 
cells with dimensions of 18 a.u. ~ 9.5 A along the di- 
rections perpendicular to the wire axis and 15 Au atoms 
are used to compute the chemisorption energies. The ir- 
reducible Brillouin zone (BZ) is sampled with a uniform 
shifted-mesh of 7 A;-points in the direction of the wire 
axis and the P point in the perpendicular directions. To 
reduce the effects of the periodic CO replicas on the elec- 
tronic bands, the number of Au atoms in the projected 
density of states calculations is increased to 105 and 6 
fc-points are used to sample the BZ. 

The ballistic conductance is evaluated with the 
Landauer-Biittiker formula, G — e^ jh T{Ep), where 
T{Ef) is the total transmission at the Fermi energy. We 
calculate the electr on transmission using the scattering- 
based approach of IChoi and Ihml 4^ generalized to the 
USPP scheme^ and implemented in the PWcond code^^ 
The scattering region is modeled with supercells consist- 
ing of a portion of the leads plus the impurity region, as 
shown in Fig. [TJ Here we use 17 Au atoms for the atop 
geometry and 18 atoms for the bridge geometry, but we 
tested longer supercells to ensure that the CO replica in- 
teraction does not affect significantly the results. Shorter 
cells with one Au atom are used to compute the complex 
band structure of the metallic leads (leftmost portion of 
the supercell in Fig. [T^). 

According to our LDA calculations,'^^ and also to pre- 
vious GGA calculations,'*^ when the Au-Au spacing in 
the pristine monatomic chain is close to its equilibrium 
value, a d band edge pinned at the Fermi level is respon- 
sible for a slightly magnetic ground-state and also for a 
large increases in the number of conductance channels, 
leading to a theoretical conductance of 3 Go . Both these 
features have never been observed in nanocontact exper- 
iments and hence can be considered a spurious effect of 
the LDA (or GGA).'^^ A non- magnetic ground state and 



the 1 Go conductance of the Au chain can be recovered 
by including atomic tips in transport calculation s'^^''*^' — 
or considering a sufficiently large Au-Au spacingi^ It 
is also possible to remove these spurious features in the 
tipless chain by including electron correlation effects at 
the DFT-I-U level, even in a very simple formulation<^ 
Here we will just consider the plain LDA results for the 
tipless chain, focusing mainly on the qualitative features 
induced by CO on the transmission function, rather than 
on the precise ballistic conductance values. Nevertheless, 
for strained geometries, where d channels are ruled out 
from the conduction, a conductance trend with respect 
to strain can be extracted from our data and compared 
with the experimental results. 



III. GEOMETRY AND ENERGETICS 

In order to understand if the adsorption site preference 
of CO may vary during different stages in the pulling pro- 
cess of the monatomic chain, we examine the chemisorp- 
tion energy and the bonding distances at different strains 
of the Au chain. Experimentally, suspended chains can 
be at most 7-8 atoms long and the chain edges are usu- 
ally attached to the leads through tips, possibly resulting 
in local variations of strain owing to the different atomic 
coordination of the atoms at edges or in the middle of 
the chain. The geometry of short suspended chains can- 
not be reproduce exactly with a simplified model not in- 
cluding the tips. Therefore, we consider here two differ- 
ent limiting cases: in the first case, we assume that the 
stress is released by the elongation of the Au-Au bond 
where CO is adsorbed (in the bridge or in the substitu- 
tional configurations), since the interaction with CO is 
likely to soften that metallic bond; in the second case, 
we distribute the strain homogeneously across the chain, 
assuming that the Au-Au bond is not significantly weak- 
ened by the presence of the molecule. The latter situation 
might be realized if the Au chain is grown on a template 
where the Au-Au spacing in the chain is imposed by the 
underlying substratci^i^ The degree of applicability of 
this model to the real nanocontact will be addressed in 
a later worki^ 

The CO/Au chain systems have been structurally re- 
laxed by optimizing the positions of C and O along x 
for the bridge and atop geometries, along z for the sub- 
stitutional geometry, and in the xz plane for the tilted 



TABLE I. Optimized distances dAu-c and dc-o (in A, see 
Fig. [1} and chemisorption energies -Bchcm (in eV) for the 
bridge and atop adsorption geometries. 



Bridge 

dAu-C dc-O -Echom 



Atop 

dAu-C dc-O -Echom 



LDA 

GGA 



1.95 
1.99 



1.16 
1.17 



-2.4 
-1.7 



1.89 
1.96 



1.14 
1.14 



-1.0 
-0.5 



bridge geometry (cf. Fig. [T]).'^^ All Au atoms are instead 
kept aligned along z with fixed positions. We first con- 
sider a chain with all Au-Au bonds at their equilibrium 
spacing, dchain — 2.51 A within LDA and dchain = 2.61 A 
within GGA, and we compare in Table U the optimized 
bond lengths and the chemisorption energy (E'chcm) for 
the bridge and atop geometries. The sum of the ener- 
gies of an isolated CO and of an isolated Au chain at its 
equilibrium spacing is used as reference energyj^ Both 
LDA and GGA calculations predict that the bridge site 
is energetically favored, by about 1.4 and 1.2 eV, respec- 
tively. Correspondingly, the C-0 bond elongation after 
adsorption is larger by about 0.01 — 0.02 A in the bridge 
geometry, as was also found in the case of Pt chains i^i^i 
As one could expected, the adsorption energies of CO on 
Au are generally smaller than those on Pt,^'^ but the dif- 
ference is much larger for the atop geometry. As a result, 
the energy difference between the bridge and the atop 
site is more marked in Au chains than in Pt chains. 

In the first limiting case of strain, we repeat the re- 
laxation of the atomic positions for different values of 
the Au-Au distance dAu-Au between the two Au atoms 
in contact with the molecule. In Fig. [21 we report as a 
function of dAu-Au the total energies of the bridge, sub- 
stitutional, and tilted bridge geometries. Three ranges of 
Au-Au bond distances can be identified: (i) when c^au-Au 
is close or slightly larger than the equilibrium spacing of 
the isolated chain, only the atop and bridge configura- 
tions are possible; (ii) when c^au-Au is very large, corre- 
sponding a Au-Au bond close to or beyond the rupture 
point, a substitutional CO in a linear configuration is 
favored;^^ (iii) at intermediate values of dAu-Au, a tilted 
bridge configuration is the most stable one. 

We note that the bridge geometry of CO is the low- 
est lying configuration not only when the Au-Au bond 
below the molecule is unstretched (^au-Au — enchain), but 
also for a wide range of ^au-Au up to about 4.2 A (cf. 
Fig. [5]). The energy minimum of the bridge configura- 
tion is at a distance c^au-Au larger than dchain, as also 
happens in Pt chainsi^ For instance, within GGA the 
minimum is c^au-Au = 2.87 A, corresponding to an en- 
ergy gain of about 0.15 eV with respect to the uniformly 
spaced configuration (dAu-Au = dchain)- The regime 
with a very elongated Au-Au bond, corresponding to 
^Au-Au > 5.5 A within GGA (dAu-Au > 5.0 A within 
LDA, see Fig. [5]), is characterized by an energy mini- 
mum for the substitutional configuration of CO, which 
is energetically favored at these very large Au-Au spac- 
ings. Within GGA, the substitutional minimum is at 
dAu-Au — 6.0 A, a distance much larger than in the LDA 
(dAu-Au = 5.16 A), and the energy as a function of dAu-Au 
is much more shallow around the minimum with respect 
to the LDA curve. Using as a reference the substitutional 
energy at dAu-Au -^ oo, the depth of the energy mini- 
mum calculated within GGA is less than 0.1 eV, while 
in the LDA it is about 0.3 eV. Given that the substitu- 
tional minimum is much deeper in Pt than in Au (more 
than three times deeper according to GGA, being about 



0.3 eV deep in Pt^), one might expect that, once CO 
is in the substitutional position, the Au-CO-Au junc- 
tion could break more easily than the Pt-CO-Pt junc- 
tion. Finally, in the short range of intermediate distances 
corresponding to 4.3 A < dAu-Au < 5.4 A within GGA 
(4.1 A < dAu-Au < 5.0 A within LDA), a tilted bridge 
configuration ("Tilted" in Fig. [5]) is preferred to the up- 
right bridge and substitutional positions. In this region, 
an energy maximum is found at about dAu-Au = 4.6 A 
(dAu-Au = 4.4 A) between the minimum of the bridge 
configuration and that of the substitutional one. 

To address the case where the strain is homogeneously 
distributed along the chain, we optimize the C and O 
position along x for a representative set of values of 
the chain spacing dchain between 2.5 A and 2.9 A. The 
chemisorption energies are obtained using as reference en- 
ergy the total energy of the isolated Au chain at the cor- 
responding inter-atomic spacing. Since it has been shown 
that a substitutional CO can only be hosted between two 
very largely spaced Au atoms, here we will consider CO 
in the bridge and atop positions only. In Table [III we 
report the optimized distances and the chemisoption en- 
ergy -Bchom- The first and second row correspond to Au- 
Au distances around the equilibrium lattice spacings of 
a pristine chain obtained from LDA (dchain = 2.50 A) 
and GGA (dchain — 2.60 A), respectively, while the third 
(dchain = 2.70 A) and the fourth (dchain = 2.90 A) de- 
scribe a chain subject to a moderate tension or near rup- 
ture, respectively. The chemisorption energies in the ta- 
ble show clearly that the bridge site is preferred also when 
the Au chain is uniformly stretched. An increase of strain 
results in larger binding energies for both bridge and atop 
geometries, but the preference towards the bridge site is 
reinforced when going to larger values of dchain- In the 
bridge geometry, the increase of i?chcm with strain can 
be ascribed to the adjustment of the Au-C-Au angle to- 
wards the optimal value. When the chain is strained 
the CO moves towards the wire axis, as seen from the 
decrease of the distance between the C atom and the 
wire axis (dc-axis), while the Au-C distance (dAu-c) stays 
almost constant. As a consequence, the bond angle in- 
creases towards the optimal value, which corresponds to 
a Au-Au distance larger than the equilibrium dchain (see 
the position of the bridge energy minimum in Fig. [21 at 
dAu-Au — 2.8 A in LDA and at dAu-Au — 2.9 A in GGA). 

Since the chain is kept straight, we find a residual 
force acting on the Au atom below the atop CO after 
relaxing the C and O positions. This force indicates a 
pronounced tendency for that Au atom to displace to- 
wards the C atom when the Au strain is not too high. 
Indeed, by relaxing also the x coordinate of that spe- 
cific Au atom we find (within GGA) that it moves away 
from the wire axis by 0.95 A when dchain = 2.60 A, but 
only by less than 0.27 A in a strained configuration with 
dchain = 2.90 A (in fair agreement with Ref. |27I ). At 
medium/high strains, ii^chem is almost stationary between 
-0.65 eV and -0.64 eV for 2.70 ^ dchain «? 2.90 A, and 
thus is only slightly larger than for the straigh chain 
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FIG. 2. (Color online) Energy of the CO/Au chain system with C and O positions optimized as a function of the Au-Au 
distance cIau-Au (see Fig. [1]). The bridge (filled circles), substitutional (filled diamonds), and tilted bridge (empty diamonds) 
geometries have been computed for several values of cZau-Au (lines through the points have been added to guide the eye), while 
the atop geometry (horizontal thick line) only at c^au-Au = dchain (see long arrow). The short arrows indicate the position of 
the energy minimum of the substitutional configuration. The reference energy (dot-dashed line) is the sum of the energies of 
an isolated CO and an isolated Au chain with c^au-Au = rfchain- This energy has also been calculated also for some values of 
dAu-Au > dchain (empty squares) . 

TABLE II. Relaxed distances (in A) and chemisorption energies (in eV) of CO adsorbed at the bridge and atop sites on the 
infinite Au chain as obtained within LDA (CCA values in parenthesis) for selected values of the inter-atomic spacing in the 
chain (dchain)- 



Bridge 



'^cha 



dc-, 



dAu-C 



dc-o 



-t-'chen 



rfAu-C 



Atop 
dc-o 



-Bchc 



2.50 1.50(1.55) 1.95(1.99) 1.16(1.16) -2.42(-1.53) 

2.60 1.45(1.51) 1.95(1.99) 1.16(1.17) -2.64(-1.72) 

2.70 1.41(1.47) 1.95(2.00) 1.16(1.17) -2.78(-1.83) 

2.90 1.31(1.38) 1.95(2.00) 1.16(1.17) -2.91(-1.93) 



1.89(1.96) 1.14(1.14) -0.99(-0.36) 

1.89(1.96) 1.14(1.14) -1.14(-0.50) 

1.90(1.96) 1.14(1.14) -1.22(-0.58) 

1.90(1.97) 1.13(1.14) -1.26(-0.62) 



(see GGA data in Table HI]). Close to the equilibrium 
Au-Au spacing, instead, we find i?chcm = —0.72 eV at 
rfchain = 2.60 A, which is about 0.22 eV larger than for 
the straight chain. Taking into account this additional 
structural relaxation, the variation of the atop Schcm 
with strain is actually the opposite with respect to that 
of the bridge i^chom. We can argue that at lower strains 
the chemisorption energy is larger because the structure 
is more close to the optimized geometry of the neutral 
AU3CO cluster, where the gold atoms form a nearly equi- 
lateral triangle with Au-Au distances of about 2.74 A. ^^ 

At this point it is interesting to compare our GGA 
data for CO adsorbed on one-dimensional chains with re- 
cent DFT calculations of CO adsorbed on Au surfaces.— 
Although the surface calculations in Ref. [5^ are at a 
much higher CO coverage and adopt a slightly differ- 
ent exchange-correlation functional (PW91 instead of 
PBE), the following conclusion can be drawn from the 
comparison:^^^ (i) distances and chemisorption energies 
for CO atop of step edge atoms in the Au(211) sur- 



face are not much different from those in the monatomic 
chain (surface: dAu-c — 1.99 A, i^chom = — 0.54eV; chain: 
dAu-c = 1.97 A, £:chcm = -0.62 eV); (ii) for CO adsorbed 
at the bridge site the differences are larger, especially for 
-Echcm (surface: dc-surf = 1.46 A, i^chem = -0.65; chain: 
dc-axis = 1.38 A, i^chcm = -1.93) (ill) larger chemisorp- 
tion energies in the chain are attributable to different 
coordination numbers of Au atoms, much smaller in the 
chain, since the comparative study of the different Au 
surfaces shows that a decrease in the local coordina- 
tion of surface atoms leads to higher binding energies 
of COj^i^i^ (iv) since the i^chem difference with respect 
to the surface case is much larger for CO at the bridge 
than for CO atop, the preference for the bridge site in 
the monatomic chain is much more marked than in the 
stepped surface (where the two sites are separated by just 
0.11 eV); (v) the C-0 bond elongation in the surface case 
(dc_o = 1.17 A for CO atop and dc-o = 1-15 A for CO at 
the bridge) is comparable to the values obtained for the 
one-dimensional chain. 




FIG. 3. (Color online) PDOS for the LDA bridge configura- 
tion. In the upper panel, we project onto the even atomic 
orbitals of C (solid red lines) and O (dashed blue lines). In 
the lower panels, the projections are on the even orbitals of 
one of the two Au atoms close to CO (solid lines). The PDOS 
for the isolated Au chain are also shown (dashed lines). 



IV. ELECTRONIC STRUCTURE 

In this section we will analyze the projected density of 
states (PDOS) of the bridge and atop geometries at the 
Au equilibrium spacing (Table |l|. We consider the pro- 
jections onto the atomic orbitals of C, O and of the Au 
atom(s) in contact with the molecule. We report only 
the projections onto atomic orbitals which are even with 
respect to the mirror symmetry through the plane con- 
taining both the CO and the chain<^ Indeed, it was al- 
ready shown for Pt chains that the even PDOS includes 
the interaction of the metal with both the a and even 
TT molecular levels of CO, while the odd PDOS involves 
only the odd tt interactioui^ There are actually slight 
differences between the even and the odd tt interactions, 
which however do not change qualitatively the conclu- 
sions obtained by examining the even PDOS only<^ 



A. CO at the bridge site 

The even PDOS for the bridge geometry are shown in 
Fig. [HI the projections on the even atomic orbitals of C 
and O (s, Px-, andp^) being in the upper part of the figure 



and those on the even atomic orbitals of one of the two 
(equivalent) Au atoms below CO (s, d322_r2, dxz-, and 
dx2^y2) in the lower part. In the PDOS projected onto 
the s and Px orbitals of C and O we can distinguish the 
peaks corresponding to the original a orbitals of CO. The 
3cr level is very low in energy (at — 23.5eV, not shown) 
and does not take part in the interaction, while the 4cr 
and 5(7 levels at — 11.2eV and — 7.7eV, respectively, have 
a moderate hybridization with Au states, as witnessed by 
the small peaks at the same energies in the Au PDOS. 
The interaction between CO and the Au chain produces 
new cr-derives states, absent in the isolated molecule, sim- 
ilarly to what was observed in a previous study of CO on 
Pt chains.— In the s PDOS of C there are three small 
additional peaks below Ep and also a very broad feature, 
centered at about 1.5 eV above Ep, which we will call 
SfJa. The latter corresponds to a cr-charge donation from 
the molecule to the metal.—"— The Sca state is mainly 
coupled with the metal s states, hence we might expected 
that it will strongly perturb the s channel transmission 
at the corresponding energies, analogously to what ob- 
served for Pt chains.— 

Among the 7r-derived states, visible in the Pz PDOS of 
C and O, the Itt and 27r* molecular states, centered at 
— 7.4eV below Ep and at 2.9 eV above Ep, respectively, 
have a small hybridization with Au states and correspond 
to small peaks in the metal PDOS. In the middle of the 
Itt and 27ra peaks, we find three separate peaks between 
—4 eV and —2 eV with distinct heights and widths. These 
states, which we call 27r^ states, arise from the backdo- 
nation of Au electrons which partially fill the formerly 
empty 27r* levels of the molecule.—"— The three sepa- 
rate contributions to the 27r^ DOS are due to hybridiza- 
tions of CO TT states with different even orbitals of Au, 
as visible from the relative intensities of the peaks in the 
Au PDOS. At larger binding energies, the two narrow 
peaks at about — 4.0eV and — 3.2eV can be associated 
to states having mostly dxz character or mostly dx^-yi 
character, respectively, while the bump at smaller ener- 
gies (about — 2.1eV) corresponds to states having pre- 
dominant s character. Thus, a donation/backdonation 
process takes place between CO and the Au chain, sim- 
ilarly to what found for Pt chainsi^ In both Pt and 
Au chains, this process leads to the formation of bond- 
ing/antibonding pairs, 5(Tb/5(Ta and 27rJ/27r*, with the 
ScTa lyh^g mainly above Ep and the 27r^ below Ep. 



B. CO at the atop site 

In Fig.m we report the even PDOS for the atop geom- 
etry, including the projections onto the atomic orbitals of 
C, O, and of the Au atom below the CO. As seen in the 
PDOS for the brigde geometry, in the s and px PDOS 
of C and O we find the Aa and Scb peaks below Ep, as 
well as the additional 5(Ta peak arising from the charge 
donation process above Ep. The 5cra peak still matches 
with a peak in the s PDOS of the metal as in the bridge 




FIG. 4. PDOS for the LDA atop configuration. The plots are 
organized in the same way as in Fig. [3] and show the PDOS of 
C and O (upper panels) and of the Au atom below CO (lower 
panels). 



geometry, but it is less broadened by the interaction and 
it is found much closer to Ep, at about 0.2 eV. We have 
already shown for Pt that the SfXa peak is correlated to 
a dip in the s-electrons transmission . ^^i^^ This dip did 
not influence the ballistic conductance of the Pt chain 
because it was located at 1.2 eV and at 2.1 eV above E-p 
in the atop"^^ and in the bridge geometry, '^'^ respectively, 
and also because of the large contribution of d electrons 
to the conductance. At variance with Pt, only s-type 
electrons contribute to the Au conductance,"^^ therefore 
one can expect that in the atop configuration the Sca 
transmission dip will also modify the ballistic conduc- 
tance of the Au chain. Finally, we notice an additional 
peak at — 3.5eV which strongly perturbs the s, d^^2_^2, 
and d^2_y2 PDOS of Au around those energies, but has 
a much smaller intensity in the s and Px PDOS of C and 
O. 

Among TT-derived even states, the distribution of peaks 
in the Pz PDOS of C and O is characterized by three 
peaks: the narrow Itt and 27r^ peaks below Ep at —7.6 eV 
and at —5.1 eV, respectively, and the broader 27r* feature 
above Ep, centered at about 2.6 eV. A single even 2tt*^ 
peak was also found at similar binding energies in the 
PDOS of the Pt chain with CO in the same position.i^i^ 
This consideration suggests that some details of the even 
TT interaction depend mainly on the adsorption geome- 



try, rather than on the metal. Indeed, in the atop ge- 
ometry even tt states are coupled prevalently with the 
dxz orbitals of the metal atom, essentially for geomet- 
rical reasons and orbital overlaps, resulting in a single 
peak. In the bridge geometry instead, there is a strong 
hybridization also with the other metal orbitals, giving 
rise to three distinct contributions (which in Pt merge 
together to give a broad peak), each characterized by a 
different strength of the coupling to d^z, dx2^y2, and s 
orbitals. 

As seen for Pt chains, the onset of these new local- 
ized states brought by the interaction with CO perturbs 
strongly the overall density distribution in the PDOS of 
the metal atoms in contact with the molecule, especially 
for states associated to narrow electronic bands, such as 
the dx2_y2 and the dxy The width of \m\ = 1 and \m\ = 2 
bands in Au chains is visibly reduced with respect to Pt 
chains and the depletion of density due to CO adsorption 
is consequently enhanced in Au. Notably, the \m\ =2 
bandwidth in Au is reduced by about a half with respect 
to Pt and the dx2-y2 and dxy PDOS suffer a much larger 
reduction of the density at energies between —1.2 eV and 
— O.leV. Therefore, at electron scattering energies be- 
tween —3.3 eV and Ep we could expect a greater reduc- 
tion in the transmission of Au chains (examined in the 
next section) with respect to that of Pt chainsi^i^ 



V. BALLISTIC TRANSMISSION 

In this section, we report the LDA electron transmis- 
sion of the bridge and atop geometries and discuss the 
effects of the uniform Au strain. As explained in Sec. |lll 
at low strains the tipless ballistic conductance results are 
affected by the spurious contribution from the d states 
at the Fermi level, hence only a qualitative discussion of 
the transmission will be presented for Au-Au spacings 
close to the equilibrium value. The conductance trend 
with respect to strain can be anyway inferred from these 
tipless transmissions because at larger strains the d chan- 
nels move to lower energies and do not contribute to the 
conduction. 

In Fig. [5l we report the transmission (solid lines) and 
the number of available channels (dashed lines) as a func- 
tion of the electron scattering energy for the bridge and 
atop geometries. We first consider the unstrained chain 
with dchain = 2.5 A (panels at the top of the figure) and 
make a connection between the electron transmission and 
the electronic structure presented in the previous section. 
Starting from the lowest scattering energies, we see that 
only one spin-degenerate to = channel is available be- 
low — 3.3eV. This channel has mostly d^z^-r'^ character 
and is well transmitted in both the bridge and the atop 
geometry, except for the smoothing at the bottom of the 
TO = band and a narrow dip at higher energies. This 
transmission dip matches with a specific peak in the s, 
d^z'^_r2, and dx2_y2 PDOS of Au, at — 3.5eV in the atop 
geometry (Fig.S]) and at — 4.0eV in the bridge geometry. 
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FIG. 5. (Color online) Total LDA transmission (solid red lines) and number of channels (dashed blues lines) for a tipless Au 
monatomic chain with CO adsorbed at the bridge or at the atop site (left and right side, respectively). Different inter-atomic 
Au spacings (dchain) are considered and the corresponding optimized distances given in Table lU are used. The arrows point 
the position of the Sua transmission dip (see text). 



in correspondence of the lower 27r^ peak (Fig. [3]). 

For energies below —3.3 eV, the transmission is slightly 
larger when CO is at the bridge and in fact, excepting 
the lower band edge region, the d^z'^-r'^ PDOS in the 
bridge geometry is less perturbed than in the atop ge- 
ometry. Conversely, at higher energies, the bridge trans- 
mission is lower than the atop transmission in the whole 
energy range of the \m\ = 1 band (—3.3 < E < 0.0 eV), 
especially for energies around — 2.2eV. In the bridge ge- 
ometry, the wide depression of the transmission centered 
at that energy can be associated with the broader 27r^ 
peak in the pz PDOS. This peak is more coupled to s 
and d^z2_j.2 states and is absent when CO is atop be- 
cause of the different orbital overlaps discussed in Sec. lIVI 
The \m\ = 2 scattering states are available only between 
— 1.7eV and — 0.9eV and are completely blocked in pres- 
ence of CO for both geometries. The \m\ =2 PDOS 
of the Au atom(s) in contact with the molecule are ac- 
tually strongly perturbed and completely detuned from 
the density of the narrow \m\ — 2 band of the chain (see, 
for instance, the d^2_y2 PDOS of Au in Fig. [3]). 

Above the Fermi level, for energies up to 4 eV only one 
spin-degenerate to = channel (mostly of s-type) is avail- 
able for transmission. In this energy region the electron 
transmission is characterized by the presence of two wide 
dips. They are present in both bridge and atop geome- 
tries but are centered at different energies which match 
exactly with the positions of 5cra and 27r* peaks in the 



corresponding PDOS. These two Fano-like features in the 
transmission can be interpreted as a consequence of a de- 
structive interference between two electron propagation 
paths: one going directly along the Au chain and another 
one passing through the corresponding resonance (5CTa or 
27r*) brought by the binding of CO to the Au chain. We 
notice that in the atop geometry the Sca transmission 
dip, pointed by the arrows in Fig. [SJ falls just 0.2 eV 
above the Fermi level, while in the bridge geometry it is 
located at about 1.6 eV. The different position of the 5(Ta 
dip turns out to be critical for the ballistic conductance 
when considering a strained chain geometry (see later) 
or when correcting for the self-interaction error of the d 
states,— since the position of the d band edge respon- 
sible for the additional (spurious) contributions to the 
conductance is very sensitive to these factors. 

The effects of strain on the electron transmission are 
described in the central and bottom panels of Fig. [SI 
with the strain increasing when going from the central 
to the bottom panels. The number of available channels 
at a given scattering energy (dashed lines) can vary with 
strain, reflecting the changes in the electronic band struc- 
ture of the isolated infinite Au chain upon stretching: (i) 
the bandwidths are progressively reduced as the Au-Au 
spacing increases; and (ii) the center of each band moves 
in energy with respect to the Fermi level. As a conse- 
quence, the positions of the d bands edges are shifted 
with respect to Ep and so are the ranges of availability 



of the corresponding transmission channels. In particu- 
lar, the number of open conductance channels changes 
as a function of the bond lengths^ at moderate strains 
(rfchain ~ 2.70 A, see central panels of the figure), the 
|?Ti| = 1 d-channels have already moved below Ep, leaving 
only one spin-degenerate s-channcl available for conduc- 
tance. However, a further increase of the strain will not 
change the number of conductance channels anymore. 

Also the transmission in presence of CO undergoes im- 
portant modifications when the chain is stretched (solid 
lines in Fig. [5]). In the bridge geometry, a moderate 
strain introduces many dips in the transmission of rn = 
states below — 2.3eV and enhances the depression due 
to 2tt^ states, which is still centered at about — 2.1eV. 
Around the Fermi energy, the d band edge resonance is 
removed and the transmission looks smoother. Above 
E-p, the ScTa dip is shifted higher in energy, following the 
shift of the corresponding peak in the PDOS (not shown 
here), and merges with the 2tt* dip to form a wide re- 
gion where the transmission is strongly suppressed. At 
higher strains (dchain = 2.90 A) the transmission at en- 
ergies below — 0.6eV degrades further, while above that 
energy it actually increases so that the ballistic conduc- 
tance slightly increases too, from 0.81 Go to 0.88 Go. In- 
terestingly, the merging of the Sca and 2tt* dips in the 
transmission leads to the partial restoring of the ideal 
transmission of the Au chain around 2.8 eV, which can 
be thought of as a result of mutual cancellation of two 
destructive interferences. In the atop geometry (right 
hand side), an increase of strain produces a progressively 
stronger transmission reduction close to the Fermi level 
and also for lower energies, down to — 2eV. At vari- 
ance with the bridge geometry, above Ep the presence 
of the SfTa and 27ra dips persists separately also at high 
strains. In particular, the Sca dip gradually approaches 
Ep, causing a progressive reduction of the conductance 
with increasing strain. Therefore, the bridge and atop 
geometries display an opposite behaviour of the conduc- 
tance with strain. 



VI. CONCLUSIONS 

In summary, we have investigated the adsorption of 
CO on a Au monatomic chain within density functional 
theory and we have studied how the electronic structure 
and the ballistic transport properties of the Au chain 
are modified by the interaction with CO. Using an infi- 
nite chain geometry we have described the energetics of 
CO in a few possible adsorption sites (bridge, atop, and 
substitutional) and, for two of them (bridge and atop), 
we have also studied its dependence on the strain of the 
Au chain. We find that the upright bridge position of 
CO is energetically favored, while a linear substitutional 
configuration becomes convenient only in presence of a 
very elongated Au-Au bond, well beyond the break dis- 
tance for the pristine monatomic chain. Thus, the site 
preference ordering for CO on Au chains is the same as 



for Pt chains r^i^i but chemisorption energies in Pt are 
substantially larger than in Au (at least 0.6 eV). The 
chemisorption energies increase with the Au strain and 
the site preference prediction towards the bridge site is 
reinforced at higher strains. 

The electronic structure of the system can be explained 
by the Blyholder model^^"— in terms of the donation 
and backdonation of electrons between the molecule and 
the metal. The adsorption process produces a broaden- 
ing and a strong hybridization of the highest occupied 
and lowest unoccupied levels of the molecule, resulting 
in the formation of bonding-antibonding pairs of 5a and 
27r'^ states: the presence of antibonding Ba^ (bonding 
27r^) states above Ep (below Ep) can be associated to 
the electron donation (backdonation). These new states 
originated from the CO-metal interaction appear as dis- 
tinct peaks in the density of states of the bridge and atop 
geometries and their binding energies, intensity, broad- 
ening, as well as the details of the coupling to the metal 
states depend strongly on the adsorption site of CO. 

In general, the electronic structure influences the trans- 
mission across the chain at all scattering energies inter- 
ested by the interaction with CO, as already seen for Pt 
chains. ■^'^■■^^ Among Au d states, those showing a stronger 
coupling with the CO levels are more poorly transmit- 
ted than the others, while s states are generally well 
transmitted except for some particular energies in corre- 
spondence of the hybridization peaks in the PDOS. Some 
dips in the transmission can be easily associated to peaks 
appearing in the PDOS at the same energies, similarly 
to the observed correspondence between the dips in the 
transmission across a molecule with a side chain and the 
molecular orbital energies of the isolated side chain. ^" In 
particular, the ba^ peak can be associated to a dip or 
a depression in the s-electron transmission. In the atop 
geometry, the ba^, peak falls very close to the Fermi level 
(E-p), while in the bridge geometry it is found well above 
Ep, hence the conductance reduction due to CO displays 
a marked dependency on the adsorption site. For Au-Au 
spacings well above the equilibrium value of the pristine 
chain, the conductance is considerably smaller in the atop 
geometry because the s channels are the only available 
for electron transport at Ep. At Au-Au spacings close 
to the equilibrium, instead, the conductance values from 
plain LDA calculations are subject to a large error due 
to the wrong position of d states^^ and do not allow a 
clear comparison between the two geometries. 

We have also studied how the transport properties of 
the chain in presence of CO depend on the strain. The 
transmission of d channels degrades at increasing strains, 
while that of the s channel is characterized by the pres- 
ence of dips which can appear, disappear or move in en- 
ergy, depending on the adsorption site. Excluding the 
low strain conductance value affected by the spurious d 
channels, ^^ the bridge and atop conductances behave dif- 
ferently as the chain gets stretched: the former slightly 
increases with strain, while the latter decreases because 
of the 5cra dip approaching to Ep. Thus, the bridge ge- 
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ometry, besides being energetically favoured, has a strain 
dependence of the conductance compatible with the con- 
ductance traces of An nanocontacts in presence of CO, 
which show a conductance increase while the Au tips are 
pulled apart. ^^ The conductance of the substitutional ge- 
ometry at the equilibrium Au spacing is found to be much 
smaller than that of the bridge geometry;^ and therefore 
it is unlikely that a tilting of the molecule from the bridge 
to the substitutional position could result in a conduc- 
tance increase upon contact stretching. 
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